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Abstract 

We explore a variant on the Mt2 kinematic variable which enables dark matter 
mass measurements for simple, one stage, cascade decays. This will prove useful for 
constraining a subset of supcrsymmctric processes, or a class of leptophilic dark matter 
models at the LHC. We investigate the statistical reach of these measurements and 
discuss which sources of error have the largest effects. For example, we find that using 
only single stage cascade decays with initial state radiation, a measurement of a 150 
GeV dark matter candidate can be made to 0(10%) for a parent mass of 300 GeV 
with a production cross section of 100 fb and 100 fb~^ of integrated luminosity. 



1 Introduction 



The Large Hadron Collider (LHC) is operational, and taking data, with expected gains in 
energy and luminosity over the next few years. One important mission for the LHC will 
be to create dark matter (DM) which appears as missing energy in the reconstructed event. 
Following a significant missing energy observation, the challenge will be to measure the 
properties of the DM candidate with sufficient accuracy to compare against cosmological 
and astrophysical constraints, such as the observed DM relic abundance and direct and 
indirect detection experiments. Thus, determining the mass of the DM particle will have 
tremendous ramifications for astrophysics and cosmology. 

Making DM mass measurements at the LHC, for example in models of supersymmetry 
(SUSY) or Universal Extra Dimensions (UED), is a difficult problem, since the DM particle 
is typically produced in pairs as products of complicated decay chains of parent particles. 
In fact, the number of states participating in the event can vary dramatically depending 
on the specific model. The identities, couplings, and masses of the particles involved in 
these processes may be unknown. Let n be the number of steps in the cascade between the 
production of the parent and the appearance of the DM child in the event. For n > 1, if all 
visible particles in the decay are detected, all masses of the parent, intermediate and visible 
and invisible child particles can, in principle, be determined uniquely (see for example pQ for 
a discussion). The simplest case of n = 1 proves to be more challenging. In Fig. [l]we show a 
schematic of an n = 1 process. We have also included the possibility that additional visible 
states are produced before the parents, which we refer to as Up-Stream Radiation (USR). In 
Sec. |2] below, we will discuss the relevance of USR for DM mass determination. Refs. [21 E] 
also study n = 1 decay chains. 




Figure 1: Schematic representation of the n = 1 class of processes considered in this work, 
with additional Up-Stream Radiation (USR). The parent particle is the state which decays 
to the visible particles and the child DM particles. 

The motivation for studying DM mass determination in n = 1 processes is many fold 
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- we mention two here. First, within SUSY or UED, n = 1 processes with additional 
USR can be important. For example, decays — )■ £^ with initial state radiation, and 
Q ^ j j i^i^ (for a sneutrino lightest SUSY particle), are of the type shown in Fig. [l| 

where £^ is a slepton, i"^ is a lepton, x° is a neutralino, g is a squark, j is a jet, is a chargino 
and z/ is a sneutrino. Although higher n chains may also be present in many models, the 
combinatoric backgrounds can make mass extraction in such decay chains complicated. By 
contrast n = 1 events are clean, and involve only two visible objects plus missing transverse 
momentum (hereafter referred to as missing energy). Also, since one will potentially observe 
n = 1 chains if one of these theories is correct, it will be useful to extract as much information 
as possible from these signals. Second, the observations of astrophysical anomalies, e.g. 
PAMELA |1] and Fermi [5], have led many to conjecture that the DM is leptophilic. Models 
which generate such signals can, for example, be constructed by connecting the DM to the 
lepton asymmetry [6], or by positing that mixed sneutrinos constitute the DM The 
simplest such dark sectors involve only a new mediator state and the leptophilic DM state, 
so that the DM is produced at a collider through the leptonic decay of the mediator. Hence, 
the study of these processes is well motivated. The reader is referred to Appendix |A] for 
more detail on models where n = 1 decay chains with USR are important. 

As shown in Appendix [B| the phase space for n = 1 processes without USR depends on 
the combination ^ = {rn^ — m^) / (2 nip) and weakly on s/(4 m^) where is the parent mass, 
nic is the child mass and a/S is the partonic center-of-mass energy. Hence, extracting /i is 
simple, while measuring proves to be more challenging. Current experimental methods for 
mass determination in events with missing energy rely on matrix element techniques. Here, 
one begins by assuming a model which implies a matrix element with additional dependence 
on TTip. Then by fitting measured differential distributions, one can extract, in addition to 
the combination /x, the overall mass scale mp by observing how quickly the event rate falls 
off with V^. 

In this paper we explore a different technique where the overall mass scale is determined 
from the transverse boosts given to the parent particles by USR. Since the boost depends 
only on mp, i.e. it is independent of the matrix element, the result is a model independent 
method for determining the overall mass scale. We explore a particular Mt2 variant proposed 
in [8], which utilizes events with USR to separately extract the parent and child masses. 
We carry out the first full scale simulation of these Mt2 based variants for dark matter 
mass determination, including detector effects, emphasizing the size of statistical errors and 
discussing various difficulties this method presents. 

The outline of this paper is as follows. We begin with a discussion of the Mt2 variable, 
and the possibility of extracting parent and child mass separately in n = 1 events with 
USR. Next we turn to a numerical analysis of this Mt2 based method and its efficiency in 
DM mass determination for a given number of n = 1 events at the LHC We then discuss 
additional sources of error beyond those explicitly contained in the previous section. Finally, 
we conclude. In Appendix A we outline some example models where this method would be 
relevant and in Appendix B we show how the phase space for n = 1 processes depends on 
the Mt2 endpoint. 
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2 Mt2 Preliminaries 



We begin by reviewing the Mt2 variable [H]. Since the LHC is a hadron coUider, the initial 
parton longitudinal momenta are unknown. Hence, only the total transverse momentum is 
constrained to be zero, and thus it becomes necessary to use transverse variables, such as 
Mt2-, a generalization of the transverse mass (see also [HI 112]). For the class of processes 
studied here (see Fig. [T]), there will be two missing particles in each event, so that the 4- 
momenta of the invisible child particles cannot be determined. Thus, only the total transverse 
missing momentum, pj^^^^, can be measured. In addition, the child particle mass, rric, is not 
known. However, a trial DM mass can be guessed, rhc, and Mj'2 formed for each event as 



max iM^^\M^^ 



(2.1) 



where the minimization is performed over trial missing momenta for the two child particles, 
^^miss(i)^ ^^miss(2)^ subject to the constraint that their sum be the total missing pt: 

^^miss(l) ^ ^^miss(2) _ ^^n,iss _ _ -'is(l) _ ^^vis(2) ^ ^3.2) 

where p^^^^^^ is the transverse momentum of the i^^ visible particle, and we are neglecting 



here the possibility of additional USR. In Eq. (2.1 ), M^^ is the transverse mass of the visible 



and child particles using the guessed missing momentum, pr^^^'^^^\ and child trial mass rric- 



M{i) If ~ 9 , { 77ivis(i) rimiss(i) ^ vis(j) ^miss(i)\ /o o\ 



where m^l is the mass of the i*^ visible particle. The energies are formed in the usual way. 



When there is no USR, there exists a value of Mt2, referred to as an endpoint, above which 
the differential cross section, da/dMT2, rapidly approaches zero, which is given by [TO] 



M^2" = f^+ V/i' + ml (2.5) 

where 

2 2 

/^^^ (2.6) 

2 TUp 

is the momentum of the invisible child in the parents' rest frame. As we show in Appendix 
[B| n = 1 chains only depend on /x up to small corrections due to the parent mass (which 
is exploited by the matrix element methods). Methods which do not capitalize on these 
corrections do not have enough information to extract both masses separately. This neglects, 
however, the potential for additional USR in the event. The USR can be in the form of jets 
coming from the initial state QCD radiation (ISR), or jets coming from the decays of heavy 
colored objects in n > 1 processes, where the decay chain ends in the n = 1 process of interest. 
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By including the USR, Eq. (2.2) no longer obtains, and instead the total momentum of the 



visible and invisible particles must be balanced against the momentum of the radiation, 



Pt 



USR 



-f miss(l) . -. miss(2) 

Pt +Pt 



Pt 



-Pt 



vis(l) 



Pt 



vis(2) 



-Pt 



(2.7) 



Now the Mt2 endpoint will depend on the upstream momentum [H [3] : 

1 1/2 



^^i-PT) + V'(M-^T)-f)' + m2^ ' . 



p2 



if rrir < rrir 



1 1/2 



p2 



if fhr > nir 



and 



2 2 

m,:, — mz 



2 m.„ 




T 



2 mp / 2 nip 



(2.8) 
(2.9) 



The functional form for the Mt2 endpoint depends on whether the test mass is larger or 
smaller than the true DM mass. Hence, there is a discontinuity in the derivative with 



respect to the trial child mass of Eq. (2.8) above and below the true DM mass, mc, giving 



rise to a kink [21 [13] in the M™2^{mc^ Pt) curve which can be utilized for extracting additional 



information beyond Eq. (2.5). In principle, given an event with a specific value for the Pt 



of the USR, one can now extract the parent and child masses. However, since one must do 
this analysis for a particular bin in Pt-, there is competition between the size of the bin - 
small bins imply small statistical samples - and the accuracy of the measurement. 

Another method was proposed in [8] , which sidesteps the problem of binning by utilizing 
the whole range of Pt- From Eqs. (2.5) and (2.8), it can be seen that M^2^ is unchanged 
by the effects of the Pt when fhc = rric- Furthermore, it has been shown [8] that 



Mr 



max / ~ 



T2 



Pt) - Mr 



max / ~ 



T2 



0) >0, 



(2.10) 



where the equality only holds when rric = rric- Thus one can construct a new variable [8] 



all events 



e {M~'''^{m,) 



Mr 



T2 



0)) 



(2.11) 



where 0(...) is the Heaviside function and M^2^^^^^'^{''^c) is the measured value of MT2{'nic)- 
It is this variable we will be minimizing to find the correct child mass. In Fig. |2| we plot 
N{ific) vs. rfic for rrip = 300 GeV and rric = 150 GeV. Since the shape is "bowl" -like, 
we refer to this construction as an Mt2 bowl. Unless otherwise specified, all events were 
simulated with the MadGraph 4.4 event generator [H], showered by PYTHIA 6.4 and 
run through the detector simulation software PGS 3.3 [16]. Note that we use the MadGraph 
default settings which defines a lepton as having pt > 10 GeV and a jet as having pt > 20 
GeV. 
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Figure 2: Example of an Mt2 bowl for 50,000 smuon pair production events with QCD 
USR. The parent mass is 300 GeV and the child mass is 150 GeV. The events were run 
through the PGS detector simulator. 

3 Mass Determination from Mt2 Bowls 

In this section we will calculate the statistical errors for child mass determination with 
bowls. Clearly, Eq. (2.9) only depends on the kinematics of the event, i.e. it is independent 



of the quantum numbers, including the spin, of the underlying particles. Then, up to small 
corrections due to the steepness of the Mt2 distribution about this endpoint, there are only 
(9(1) differences in the bowls around the minimum for different parent spins. Hence, we can 
study the effectiveness of this variable for a wide variety of models by only scanning over 
the masses of the parent and child particles. We take 



mp = 100 GeV, 
rrip = 300 GeV, 
TUp = 500 GeV, 



= 25, 50, 75 GeV, 
m, = 75, 150, 225 GeV, 



rric = 125, 250, 375 GeV 



(3.1) 



as our benchmark parameters. 

For reference we provide the overall cross section for these benchmark models in Table [T| 
where we have assumed that the production occurs via electroweak processes, including the 
effects of QCD ISR. Neglecting diagrams which involve additional new-physics states, the 
overall rates only depend on the spin of the parent up to 0{1) factors due to the choice of 
SU{2) X f/(l) representation. For reference, the scalar example process is pp i~ — )■ 

£~ x° where £^ is a slepton and x° is the lightest neutralino. This is the process we 



6 



simulate for our benchmarks with QCD ISR. For reference, a fermionic example process is 
PP XT ^ i~' vv* where is a chargino and i> is a sterile sneutrino. For some 

details of these and other models which have n = 1 processes, see Appendix |A} Our results 
below will be given in terms of the number of events before cuts, so Table [T] can be used to 
estimate the reach of actual models. 



nip 


'^scalar 


C^fcrmion 


100 GeV 
300 GeV 
500 GeV 


0.4 pb 

9 X 10-3 pb 
10-3 pb 


20 pb 
0.4 pb 

6 X 10-2 pb 



Table 1: Cross sections for electroweak pair production of parent particles with various 
masses and spins including the effects of QCD ISR. We neglect any t-channel processes 
involving additional states. 



There are also models which have more complicated decay chains but can be interpreted 
as n = 1 processes with additional USR. For example, one can have new colored objects 
which decay to jets and the parent particle. As long as the USR can be distinguished from 
the decay product of the parent, our method is applicable. This will improve the prospects 
for this method dramatically since the overall rate will increase due to colored production 
instead of electroweak production, and additionally the majority of events will have very 
hard Pt for the USR. Hence, we also choose a set of benchmark models with colored objects 
up-stream with masses 

mcoi = 600 GeV, rup = 300 GeV, = 150 GeV, 

mcoi = 1000 GeV, rUp = 300 GeV, me = 150 GeV, 
mcoi = 1400 GeV, nip = 300 GeV, = 150 GeV, (3.2) 

where rricoi is the mass of the colored state which decays to the parent particle and jets. The 
example process we will simulate for these benchmarks is pp ^ q q ^ j j x~ d-^ ^~ v v* 
where g is a squark, is a chargino and i> is a sterile sneutrino. While one might be able 
to use additional handles from viewing such events as n = 2 processes (instead of n = 1 with 
USR), here we wish to examine only the effect of harder USR from colored particle decay 
on the error for DM mass determination in n = 1. For additional models of n = 1 processes 
which can be produced in the decays of colored states see Appendix \K\ 

As discussed above, the additional radiation shifts all events, including those near the 
Mt2 endpoint. For reference, the radiation distributions for our benchmark models are 



shown in Figs. |3j and |4j From Eq. (2.8), the correction to the Mt2 endpoint due to USR 
is of the form Px/mp. Hence, the Pt distribution of jets determines how well the parent 
and child masses can be extracted separately. From Fig. [3} we see that heavier parents lead 
to harder P^ distributions due to the larger recoil occurring from production of a heavier 
state. However, since the correction to Mt2 goes as l/m^, this enhancement is tempered 
by the parent mass. In addition, heavier parents have smaller production cross sections 
(see Table [T]). Hence, assuming they can be seen above the backgrounds, lower parent mass 
states give rise to more defined bowls. The trade-off between background rejection, which 
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is optimized for high masses, and the quahty of the Mt2 bowls, which is optimized for low 
masses due to the dependence on the ISR, leads to a sweet spot in the range of C(200 GeV) 
to (9(500 GeV), with significant dependence on the spin of the parent. In the cases with 
colored states upstream this tension is alleviated since now the Pt distributions are harder 
and the production cross sections are larger, as in Fig. |4j 




100 200 300 400 500 

Pt [GeV] 



Figure 3: Pt of the hardest jet for slepton events with only QCD ISR. The blue dotted line 
is for nip = 500 GeV, the red dashed line is for nip = 300 GeV and the yellow sohd hne is 
for nip = 100 GeV. Note from Eq. (2.8) that the correction to Mt2 due to USR is of the 
form PT/nip. 



The effects of the radiation on the Mt2 endpoint are shown by plotting N{nic) as a 
function of nic in Fig. |2} for 50,000 smuon pair production events with two muons and 
missing energy, with no background events (also see Fig. 11). As we will show in the next 
section, the backgrounds can be very efficiently cut away, and will be insignificant near the 
Mt2 endpoint (see Sec. 4.2). In what follows, we will present statistical error bars on the 
DM mass determination using the Mt2 bowl and will discuss in detail various sources of 
error and their effect on this analysis. 



3.1 Statistical Analysis of Mt2 Bowls 

Contributions to adjacent bins in the Mt2 bowls from the same events imply that it is 
inappropriate to use simple a/ZV statistics in computing errors. Removing one event from a 
given bin in the distribution can in principle remove one event from each bin. Therefore, we 
utilized the well-known "bootstrapping" method to do the statistical error analysis. We 
employed the following method when doing this. We begin by generating a sample of 
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Pt [GeV] 



1500 



Figure 4: Pt of the hardest jet with new colored state dominating the USR. Specifically, 
these colored states are squarks. The blue dotted line is for rricoi = 1400 GeV, the red 
dotted line is for rricoi = 1000 GeV and the yellow solid line is for mcoi = 600 GeV. We have 
fixed nip = 300 GeV in all cases. From Eq. (2.8), the correction to Mt2 due to USR is of 
the form Px/rUp. 



C(100,000) signal events (we take y/s = 14 TeV). From those 100,000 events, we choose 
a subset of size Ng^ents, and make 100 independent random selections of Nevents events from 



the original data set. Then for each of these selections we calculate N{rfic) using Eq. (2.11). 
This gives us a random sampling of bowls for a given number of events. Since there is often a 
degeneracy of minima for each of these random bowls, especially for a low number of events, 
we take the geometric mean of these multiple minima to give us an average minimum for 



each bowl. Note that we do this assuming the theoretical value of M^2^ (see Sec. 4.3). 
Finally, we find the mean and standard deviation of these 100 average minima. To find the 
standard deviation we used the formula ^(xj — a;mean)^/(^— 1) and checked to confirm that 
this corresponds to 1-cr error for a Gaussian distribution to good approximation. 

This method allows for a statistical sampling of the distribution of possible bowls for a 
given number of events. We present our results as a function of Ai'events, the number of events 
before any cuts are made. Note that the events which contribute to the bowl have very 
special kinematics which allow them to go beyond M^2^ - the overwhelming majority of 
events will not have any bearing on the mass determination. Hence, cuts designed to remove 
backgrounds will not cut away these special events which contribute near the minimum 
of the bowl where the DM mass determination occurs. This is an expectation we 
check explicitly in the next sectionj^] Also note that by working with the mean we will 



*This assumption is not true when MJpf^ « mw as in the case of, for example, nip = 100 GeV and 
rur = 25 GeV. 
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systematically underestimate the DM mass due to the asymmetric shape of the bowl. This 
asymmetry is due to the shape of the Mt2 distribution near the endpoint as a function of 
rhc - the slope becomes steeper as rhc is taken larger. The events used for the bowls were 
generated using the PGS detector simulator so that they do include detector effects which 



also adds to the consistent underestimates. As we discuss in Sec. 4.1, detector simulations 



must be utilized to determine the required correction to account for this off-set. Further 
sources of error are discussed below in Sec. HI 

In Figs. [5] - 18} we show the statistical error bars for the DM mass determination for a 
given parent and child mass combination as a function of the number of events before cuts. 
Note that for a given child mass, the error bars grow smaller as the DM mass approaches 
the parent mass, due to the width of the minimum of the bowl. This occurs because the 
minimum of the bowl becomes more well-defined as the Mt2 distribution becomes steeper. 
The error bars grow smaller as N events grows larger, but not as quickly as l/V-^events- This is 
because events contribute to multiple bins so that errors from adjacent bins are correlated. 
Also notice that error bars in Fig. [8] are much smaller for a given iVevents than those in Fig. |6] 
for rric = 150 GeV. The error bars are also smaller for larger values of mcoi- This is due 
to the enhanced Pt of the USR as shown by comparing Figs. [3] and |4j In what remains we 
will discuss the various additional errors and will argue to what degree we expect them to 
degrade the results. 
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Figure 5: Mean and ±l-cr statistical error bars for a DM mass measurement as a function 
of the number of signal events before cuts. The only source of USR is initial state radiation. 
The process we simulated is electroweak smuon production. The error bars will improve by 
0{1) for fermionic parents. The parent mass is 100 GeV and the child masses are 75 GeV 
(green), 50 GeV (blue) and 25 GeV (red) from top to bottom. The dashed lines show the 
actual child mass. Note that detector effects have been simulated for the underlying events 
and that the DM mass measurement systematically undershoots the actual value on account 
of these effects. 
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Figure 6: Same as Fig. [5] except that the parent mass is 300 GeV and the child masses are 
225 GeV (green), 150 GeV (blue) and 75 GeV (red) from top to bottom. As explained in the 
text, we find that cuts designed to eliminate the background will not change these results. 




Figure 7: Same as Fig. [5] except that the parent mass is 500 GeV and the child masses are 
375 GeV (green), 250 GeV (blue) and 125 GeV (red) from top to bottom. As explained 
in the text, we find that cuts designed to eliminate the background will not change these 
results. 
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Figure 8: Same as Fig. [6]except that the dominant source of USR is new heavy colored states. 
These colored states are squarks which produce chargino parents and jets. The parent mass 
is 300 GeV and the child mass is 150 GeV for all three cases. The mass of the colored 
objects are 1400 GeV (green), 1000 GeV (blue) and 600 GeV (red) from top to bottom. As 
explained in the text, we find that cuts designed to eliminate the background will not change 
these results. 



4 Sources of Error 



The results of Figs. [5]-|8]only incorporate statistical and detector effects. In this section we 
argue that the errors we have included in our analysis are a realistic estimate of the precision 
with which the DM mass can be extracted from simple cascade decays. We further qualify 
the additional sources of error below. 



4.1 Detector Effects 



With the inclusion of detector effects, the events at the Mt2 endpoint become smeared 
out. This implies that some events which do not have the correct kinematics to make a 
contribution to the bowl can have Mt2 > M^f^. This leads to a degradation of the minimum 
of the bowl. Since the Mt2 distribution is steeper for larger test masses, this degradation 
will tend to contribute to a larger underestimate of the DM mass. This is the reason for the 
systematic under-shooting of the DM mass in Figs. |5] - [8j To illustrate this effect we have 
generated the analog of Fig. [6] for parton level events as shown in Fig. [9j Note that the l-cr 
error bars overlap with the actual DM mass except in the case where = 75 GeV since 



here the bowl is essentially flat below rric ~ 75 GeV (see Fig. 11 ). Hence detector simulations 



would have to correct for this systematic effect in any real DM mass measurement. 



After generating bowls using the parton level events, the value N{mc) (see Eq. (2.11)) at 
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the minimum is ~ 0. For the same bowls, but with detector effects, the value N{mc = rric) 
is no longer - for (9(100,000) events, N{fhc = rric) ~ (9(100). Hence one can attempt 
to clean up the bowl by removing the events from the data sample which contribute at the 
minimum. This will increase the steepness of the bowl and might be helpful in minimizing 
the error since all removed events are guaranteed to be pathological. However, since this 
cleaning process does not change the minimum, this will not change the error bars presented 
above. 

Parton Level Events 
mp=300 GeV 



225 




500 700 1000 1500 2250 3500 5000 7000 10000 15000 22500 

^events 



Figure 9: Same as Fig. [6] except that the underlying events are parton level. 
4.2 Background Contamination and Cuts 

In this section we will argue that a generic set of cuts designed to remove backgrounds will not 
degrade the minimum of the Mt2 bowl and hence will not affect our conclusions. Motivated 
by the choices taken in [17], we have analyzed the following cuts for illustration, which are 
relevant for di-lepton events with jets and missing energy {i.e. slepton pair production): 

1. Require 2 opposite sign, same flavor leptons (e or /i). 

2. Hardest lepton: pt > 40 GeV. 

3. Second hardest lepton: px > 30 GeV. 

4. p^''' > 100 GeV. 

5. A veto: the invariant mass of the two leptons, mu, must not lie in the range 
80 GeV < mu < 100 GeV. 

6. No 6-tagged jets. 
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While cuts should be tailored to the particular model under consideration, these are fairly- 
generic, and will serve to illustrate the point that our results are not significantly degraded by 
background removal. We also explored the effect of a cut on Mt2 by requiring Mt2{'^c = 0) > 
100 GeV. These cuts will be very efficient for eliminating standard model (SM) backgrounds, 
the worst of which is W~ plus jets, where the bosons decay leptonically. In 
particular, this di-boson process is dominated by 1 1 production. 

An Mt2 cut on the 1 1 background is a powerful discriminator, and in many cases it will 
have no effect on the DM mass determination. To see this, first note that the t i background 
falls into the same class of n = 1 processes we have been studying already, with the tops 
as the colored particles leading to hard USR, the as parents and the neutrinos as 
children. Since the child is a neutrino, rric = 0, and the minimum of the bowl will occur 
at rhc = 0. Then (neglecting detector effects which will only add a small perturbation) 



the ti background will be largely eliminated for an Mt2 cut of (9(100 GeV). In Fig. 10 
we plot this Mt2 distribution including detector effects. Clearly, there is an endpoint at 
niw- The cross section for 1 1 — )■ 6 5 /i~ /x"*" z/^ z/^ is 5 pb. Then starting with a 100,000 event 
sample, the cuts 1-6 described above reduce this background to 0.065 ± 0.002 pb. Then the 
Mr2(^c = 0) > 100 GeV cut eliminates all remaining events. In this way, the worst of the 
SM backgrounds can be easily removed for nip ^ 100 GeV. 
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Figure 10: Mt2{i^c = 0) distribution for tt events where we have treated this as an n = 1 
process where the 6-jets are USR and the are the parent particles. This plot is made 
before cuts and we have included detector effects. There is an endpoint at mw since the 
child, i.e. the neutrino, mass is zero in these events. 



In Fig. 11 we have plotted a series of Mt2 bowls before and after this set of cuts 



to check that the signal in the DM mass determination region of the Mt2 bowl is not 
degraded. For rrip = 100 GeV there is a significant degradation of the bowl. However, 
for this value of m„, there will be tremendous difficulties disentangling the signal from the 
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W~ background since they have very similar Mt2 endpoints. For the models with heavier 
parents or with additional colored states producing hard USR, the minimum is maintained 
for these cuts. Additionally, the Mt2 cut has no effect on these plots (excluding the example 
with rup = 100 GeV). We also checked that this statement is robust under variations in the 
cut parameter choices made above. 



4.3 Variation in M^^"" 



In generating Figs. |5]j8]we assumed that the endpoint has been measured precisely and 

matches the theoretical value. In [S], another Mt2 based variable, Mt2±, was introduced, 
which is the projection of along the direction perpendicular to the USR. They show 
that the endpoint of this distribution is independent of the USR momentum and identical 
to M^2^{rhc, 0) endpoint. Hence, even in cases with large USR, it is possible to extract the 
required input to construct the bowls. 

However, the level of accuracy with which M^2^{mc, 0) can be measured depends on 



detector effects. For the purposes of illustration, in Fig. 12, we show how the Mt2 bowl 
is degraded as one varies the M^2^ endpoint by ±2% and ±5% for nip = 300 GeV and 
rric = 150 GeV. For variations on the order of —5% the minimum is shifted by a non-trivial 
amount and can even disappear in some cases. For overestimates of M^2^ of order 5%, the 
width of the minimum becomes much broader than the statistical error bars presented above. 
Therefore, it is crucial to the success of this method that an accurate measurement of the 
Mt2 endpoint be made. On the other hand, the steepness of the bowl around the minimum 
is maximized for the correct choice of M^2^. By combining this observation with the direct 
measurement of the endpoint, the accuracy with which M^2^ could be determined would be 
improved. Other variables, such as that suggested in [] could be used to determine the Mt2 
endpoint. The accuracy with which this can be done is left for future work, though it can 
likely be done with high precision due to the larger amount of statistics available than for 
the bowls. 



5 Discussion and Conclusions 

In this work we studied the possibility of using n = 1 single stage cascade decays to measure 
the DM mass at the LHC. We have argued, using the particular Mt2 variant of [8J, that 
if a signal is observable and backgrounds can be eliminated, it is possible to make (9(10%) 
measurements of the DM mass with (9(10,000) events before cuts for optimal values of rUp 
and rric- We have shown that this requires a precise determination of M^2^{mc, 0). 

In [TH] the matrix element technique was used to ascertain how well the neutralino mass 
could be measured in an n = 1 squark decay for a benchmark model with a parent mass 
of 561 GeV and a child mass of 97 GeV. Using parton level events so that jet smearing 
effects, etc., are not considered, they found that with 3000 events before cuts only an upper 
limit on the child mass could be determined and with 7500 events a measurement could be 
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Figure 11: bowl for 25,000 (10,000) slepton (squark) pair production events which 

give jets, two muons and missing energy. The bowls on the left column only have QCD 
ISR. The bowls on the right have additional colored states which dominate the USR, and we 
have taken = 300 GeV and rric = 150 GeV for these cases. Note that the cuts preserve 
the minimum in all cases except rup = 100 GeV. Additionally, when one does a cut on 
Mr2(?^c = 0), the bowl will be unaffected as long as this cut is taken below M^2^ for the 
bowl in question. 
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Figure 12: Plot of Mt2 bowls allowing for variations in the endpoint of ±5% and 

±2%. All bowls are made with 50,000 smuon pair production events before cuts. For clarity 
we have not simulated detector effects for these events. 



made with an (9(100%) error bar. This can be compared with our Fig. [7|for the benchmark 
nip = 500 GeV and = 125 GeV[] We find that with 3000 events we can make an O{70%) 
determination and for 7500 events error bar goes down to 0(50%) once the correction for 



detector effects is applied as described above in Sec. |4.1[ Hence, the methods seem to be 
competitive, but ultimately a detailed study will be required to determine which will lead 
to the best DM mass determination. 

Finally, we would like to emphasize the model independence of these results, even when 
there are complicated cascade decays. A large class of events can be interpreted as n = 1 
processes with USR. All that is required is that the only missing energy in event is produced 
at the end of the chain as the result of the decay of an on-shell parent, and that the USR be 
distinguishable from the decay product of the parent. When this isolation is possible {e.g. 
the two photon plus missing energy signal of some gauge mediated SUSY breaking models) 
our results can be applied up to differences due to detector effects. 
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Appendices 



A Benchmark Models 

The n = 1 events studied here are the simplest class of events at the LHC which involve 
the DM. Perhaps the most commonly studied of such processes is pp q q j j 
such models, however, one expects there to be higher n processes present as well which will 
give additional kinematic information. In this appendix we will outline examples of n = 1 
process with scalar, fermionic and vector parents. Estimates for the electoweak LHC cross 
sections for these models are given in Table [Tj 



A.l Scalar Parents 

We begin by motivating scalar parents. Recently, a wave of leptophilic DM models have 
been proposed to explain measured cosmic ray anomalies. A non-supersymmetric example, 
which is additionally motivated by the baryon-DM coincidence, can be constructed by simply 
extending the SM by two additional fields: a new Higgs doublet, H', and a leptophilic DM 
state, X, interacting via [T9] 

C = XLH' + mxXX. (A.l) 

When the additional term 

AC = \{H^H'f + h.c. (A.2) 

is added to the Lagrangian, where H is the SM Higgs doublet, and the H' is integrated out, 
the effective operator 

_X^LHLH 

Casym — (A-3) 

is generated, where M is the effective suppression scale. This operator transfers the lepton 
asymmetry to the DM sector, so that the DM density is set by an asymmetry and not thermal 
freeze-out. Also note that such leptophilic DM candidates can be viable as an explanation 
for the observation of an excesses of cosmic ray positrons by the PAMELA experiment [6]. 
Although the DM would be asymmetric {i.e. mostly X) when its density freezes in, that 
asymmetry could be erased through Majorana mass terms for X and X. Then in the universe 
today, XX — )■ may give rise to significant cosmic ray positron signals. 

The DM would be created at the collider through the electroweak production of the H', 

pp^H'H'^XXtr. (A.4) 



However, production rates for pp H' H' — >■ X X £~ will be low (see Table |lj). While 
these events could be extracted from the large di-boson background with high luminosity, 
DM mass determination will be difficult. 

Note that this process is identical to the electroweak pair production of sleptons (see [IT] 
for a study which determines how feasible it is to find these processes at the LHC), 

pp^£+r ^£+r (A.5) 
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A. 2 Fermionic Parents 



For an example with fermionic parents, we turn to a model which is embedded within the 
MSSM. Introduce a superfield DM candidate, X, with the quantum numbers of a sterile 
neutrino. Then the active sneutrino can mix with scalar partner for X, X, leading to 
mixed sneutrino DM. In [7], X has been shown to be a viable DM candidate. At the LHC, 
electroweak production can go through 

PP^x^x' ^xx*e-^r, (a.6) 

where x^ is a chargino. Since the parent particles are fermions instead of scalars, the 
production rates are larger (see Table ig. 

With a slight modification, these classes of DM models can be related to the lepton 
asymmetry. One can add a new pair of electroweak doublet superfields, D and D, and a new 
superpotential term, 

AW = mDDD + XX DHu + ViL.DX + mxXX, (A.7) 

where mo is the mass for D, mx is the mass for X and A is a new yukawa coupling. 
Integrating out these doublet states results in the lepton number transferring operator 

where M is the effective suppression scale. This operator can be used to generate the relic 
density. The production at the collider then goes through the electroweak production of the 
fermionic D: 

pp^ D+ b- ^ xk* t r. (A.9) 

Production rates in all these fermionic parent models can further be enhanced by em- 
bedding the n = 1 process into squark decays: 

pp-^ qq^ X^xT n i'i>*i^i~ jj (A.IO) 

As described above (see Fig. |4]), this will lead to a much harder USR distribution, which in 
turn will imply better DM mass determination. 

A. 3 Vector Parents 

Lastly, we note that within UED models, pair production of vectors gives rise to similar 
signals. For example, 

pp W^^^- 1 1 z/(i)z/W, (A.ll) 

where is a KK ly-boson, z/^^^ is a KK neutrino is an n = 1 chain. This process can 

similarly be embedded in the decay of new colored states, which gives rise to harder USR: 

pp^ Q^^^Q^^^ lyW+W^W- jj ^ tr z/(i)i?« J j, (A.12) 

where ' is a KK quark. Note that if z/(i) is the DM its mass is restricted to be greater 
than C(50TeV) by direct detection experiments [20] . 
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B Phase Space Dependence on Mt2 



To show the phase space dependence on Mt2 and the overall scale mp, we will assume that the 



parents are produced on-shell so that the 2 — )■ 4 production in Fig. 13 can be approximated 
by the 2 — 2 cross section a2^2 and parent particle decay width V. 




Figure 13: Process considered in this section. The proton momenta are g^, the parent 
momenta are /cj, the visible momenta are and the child momenta are Pa- 

We begin by simplifying the general 2 — ?■ 2 differential cross-section and 1 — 2 differential 
decay width. Throughout the calculation we will drop overall constants since they do not 
contribute to the normalized distributions. The 2 — )■ 2 differential cross-section is given by 

= l-M.r (2 vr)^5^(gi + - - h) (B-l) 

4kilcMVS 2Ei (2 7r)3 2E2 (2 7r)3 

where a/^ is the parton center-of-mass (CM) energy and Ei is the energy of the i^^ parent. 
Integrating over k2 in the CM frame to eliminate 6^{qi + 92 ~ ~ ^2) gives 

da2^2 oc \ \M^^ -J— 6{^s -E2- El) d'h, (B.2) 
where ^2 = —ki. Similarly, we simplify the 1 — > 2 differential decay widths 

dr. ^ ^ \Mr\' (2M - p., - P.)^g5^^|^. (B.3) 

where c and v stand for child and visible, respectively, and i = 1,2. Integrating over pc- to 
eliminate 6^{ki — Pa — PvJ gives 

dF, (X ^ J ^ \Mrf 5{E, - E,^ - d=^p,„ (B.4) 



Ei Ec- Ey. 
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where the (5-function enforces = ki —pv,- Since for 1 — i- 2 decays the summed and squared 
matrix elements |A^rJ^ are only functions of the masses, they will not contribute to the 
normalized distributions. We drop these factors from here forward. 

Convolving the differential parent decay width with the differential 2 — > 2 cross section, 
and again dropping overall constant factors gives 

da = dcr2-^2 dFi dr2 

S J^ci ^vi ^C2 ^V2 

Define cos Pi to be the angle between the visible particle momenta and the parent particle: 

h -pv, = hPv, cos A. (B.6) 

Rewriting the phase space delta functions so that cos/3i, cos/32, and ki/nip are the 
integration variables, the integrand takes on a more revealing form 

dcr2^2dridr2 oc d (—) dp.-, dp^,, df^i d(cos dfcos/?^,,) d0., d0t,„ 

\mpj 

X '-^^^ jiOu e,,X2 ■■ Oi, (3v„P.2) s (— - Vs^ 



X ^|co.A-^i^^ U co.ft-^i^||S^|, (B.7) 

Pv-,VS - I J \ Pv2Vs - 1 

where J7'(...) is the Jacobian for converting from integration over the 6 angles to the /3 angles, 
which does not depend on any mass parameters, n is defined as in the main body of the 



paper (see Eq. (2.6)), and 



4m2 



(B.8) 



Hence, Eq. (B.7) shows that the phase space only depends on the parent and child 
mass through the two functions s and fi. Note that when integrating over the parton 
distribution functions (PDFs), the factor 1/(1 — s)^ will cause the differential cross-section 
to be dominated by values s ~ 1, which corresponds to threshold production of the parent 
particles. Therefore, for a trivial cross section matrix element, the differential cross section 
only depends on the rrip and m^. through the combination 2 fi = M^2^{rhc = 0) = (m^ — 
ml)/mp. Note that in some cases, A^^- will depend explicitly on nip, or masses of particles 
being exchanged in the corresponding Feynmann diagram, causing slight deviations in the 
normalized distributions for the same M^2^, but different parent mass. 

In Fig.jBj we plot various normalized distributions with the same M^2^(mc = 0) endpoint. 
With the exception of a weak dependence on the parent mass, due to s dependence, the 
distributions look virtually identical. Note that it is this weak s dependence that the matrix 
element methods seek to capitalize on. 
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Figure 14: Various distributions for points with the same Mt2 endpoint for pp ^ — ^ 
As shown in Appendix [b| the data only depends on M^2^ with a shght 
variation due to s/(4mp). The distributions plotted are the invariant mass of the two visible 
particles (upper left), the total missing transverse energy (upper right), the total transverse 
momentum of the visible particles (lower left), and Ai? = y^Ar7M-~A0^ between the two 
visible particles (lower right). 
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